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Summary 

^Predictions  of  the  attached  turbulent  boundary  layer  on 
three  axisymmetric  bodies  were  obtained  usina  integral  and 
finite-difference  methods.  The  integral  method  was  developed  by 
Myring  and  the  finite-difference  method  by  Wang  &  Huang.  Both 
methods  take  account  of  viscous-inviscid  interaction  effects  and 
the  finite-difference  method  also  takes  account  of  the  variation 
of  static  pressure  across  the  boundary  layer.  Comparisons  with 
experimental  data  show  both  methods  to  be  of  similar  accuracy 
over  all  but  the  last  few  per  cent  of  the  body  length.  However, 
the  finite-difference  method  is  marginally  superior  in  this 
region  and  gives  more  detailed  information,  including  the 
individual  velocity  components  and  the  Reynolds  shear  stress. 
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LIST  OF  SYMBOLS 


1  2 

Skin-friction  coefficient  =  (wall  shear  stress) /^pU^ 

H  Shape  factor  =  A*/0 

L  Body  length 

RO  Body  radius  at  axial  distance  X  from  nose 

R  Radial  distance  from  axis 

U.V  Tangential  and  normal  components  of  mean  velocity 

U' ,V'  Tangential  and  normal  turbulent  velocity  fluctuations 

Velocity  at  edge  of  boundary  layer 

UTNF  Freestream  velocity 

UX.UR  Axial  and  radial  components  of  mean  velocity 

X  Axial  distance  from  nose 

Y  Normal  distance  from  body  surface 

a  Angle  of  tangent  to  body  surface  with  respect  to  axis 

6  Boundary  layer  thickness 

S*  Boundary  layer  displacement  thickness 

*  f® 

A  Displacement  area  =  (1-U/U  )  (RO  +  Y  cosa)  dY 

Jo  ^ 

0  Momentum  area  =  U/U  (1-U/U  )  (RO  +  Y  cosa)  dY 

Jo  e  e 


INTOODUCTION 

Current,  ARE  propul  snr  design  methods  for  underwater 
vehicles  require  a  knowledge  of  the  boundary  layer  flow  in  the 
reoinn  of  the  propulnor  on  the  unpropelled  vehicle.  A  typical 
vehicle  conf iguration  often  comprises  a  nnn-axi symmetri c  hull 
shape  fitted  with  control  fins,  which  results  in  a  highly  three- 
dimensional  flow  field  with  vortex  sheddlna  and  high  levels  of 
turbulence.  Although  such  flow  fields  can  be  measured  on  scale 
models,  the  Reynolds  numbers  involved  are  usually  considerably 
lower  than  the  full-scale  value,  sometimes  by  a  factor  of  100. 
However,  axl symmetric  calculation  methods  can  be  used  to  give 
design  Information  and  estimate  the  boundary-layer  Reynolds 
number  scaling  effects. 

An  integral  method  for  predicting  the  turbulent  boundary 
layer  on  axlsymmetric  bodies,  developed  by  Myring  fll,  has  been 
in  use  at  ARE  for  several  years.  The  method  was  designed 
primarily  for  predicting  drag  and  has  been  shown  to  give  good 
results  for  a  wide  variety  of  body  shapes  U?.!.  There  is  also 
reasonable  aareement  with  experiment  of  the  bounda ry- layer 
integral  parameters  and  velocity  profiles,  which  are  required  for 
propulsor  design  purposes,  on  body  shapes  with  small  longitudinal 
curvature  and  tail  cone  semi  angles  of  less  than  about  20°.  For 
body  shapes  with  fuller  sterns,  the  agreement  is  less  favourable 
with  the  shape  parameter  in  particular  beina  under-predicted  at 
the  tail  end  of  the  body  where  the  boundary  layer  may  be  close  to 
separa  tion . 

Tlie  present  investigation  aims  to  compare  the  Myring  method 
with  a  more  sophisticated  finite-difference  method  developed  by 
Wang  A  Hviang  r?,4,51  and  described  in  Section  2  below.  Accurate 
measurements  of  the  boundary  - 1 ayer  velocity  components  and 
Reynolds  stresses  using  hot-wire  or  laser-Doppler  anemometry  are 
now  available,  and  these  are  compared  in  Section  3  with  pre¬ 
dictions  using  both  integral  and  f i ni te -di f f erence  methods  for 
three  representative  body  shapes. 


theciry 

Roth  the  integral  and  the  f ini te -d i f f erence  method  use  the 
same  technique  based  on  the  displacement  body  concept,  i.e.  the 
static  press\ire  or  velocity  at  the  edge  of  the  boundary  layer 
f required  by  the  boundary- 1 ayer  calculation)  is  the  same  as  the 
surface  pressure  or  velocity  in  an  inviscid  flow  over  the  body 
surface  to  which  the  displacement  thickness  has  been  added.  This 
moans  that  the  bounda ry- layer  and  potent ia 1  - f low  calculations  can 
be  coupled  together  in  an  iterative  cycle.  7T»e  potential  calcu¬ 
lation  in  both  methods  is  based  on  the  source  panel  method  of 
Hess  ft  Rmith  161,  The  main  difference  between  the  two  methods  is 
in  the  boundary- layer  calculation.  In  the  finite-difference 
method  the  momentum  eqviation  is  solved  directly,  whereas  in  the 
integral  method  It  is  formally  integrated  across  the  boundary 
layer  to  give  the  momentum  integral  equation,  which  is  then 
solved  empirically. 


The  Myring  integral  method  proceeds  as  follows:  the  boundary 
layer  is  calculated  using  an  initial  guess  at  the  external 
pressure  distribution,  "nie  location  of  transition  to  turbulence 
miKst  be  prescribed  explicitly  by  the  user.  The  potential  flow  is 
then  calculated  over  the  “displacement  surface"  and  the  predicted 
pressure  distribution  is  used  as  the  input  to  the  next  boundary- 
layer  calculation.  The  calculation  is  extended  into  the  wake  and 
suitable  end  conditions  are  imposed.  The  profile  drag  is 
estimated  from  the  momentum  deficit  in  the  far  wake  and  the 
skin-friction  drag  by  integrating  the  local  skin-friction 
coefficient  over  the  body  surface.  An  iterative  scheme  results 
which  stops  when  the  iterated  values  of  drag  have  converged, 
usually  after  about  six  iterations.  Profiles  of  the  tangential 
velocity  along  a  normal  to  the  hull  can  be  determined  by  assuming 
a  power-law  profile.  The  method  has  been  adapted  to  include  a 
simple  actuator-disk  representation  of  a  lightly-loaded  propeller 
C7n. 


The  Wang  S  Huang  method  is  based  on  the  finite-difference 
solution  of  the  boundary- layer  equations  developed  by  Cebeci  & 
Smith  Ca3.  The  location  of  transition  can  be  prescribed  or 
predicted  using  a  simple  empirical  method,  and  the  effects  of  a 
tripping  device  can  be  incorporated.  A  simple  eddy-vi scosity 
model  of  turbulence  is  used,  which  is  based  on  the  traditional 
mixing-length  approach.  The  flow  in  the  wake  is  calculated  usincr 
an  integral  method  proposed  by  Granville  [91.  The  displacement 
surface  is  determined  from  the  boundary- layer  and  wake  calcu¬ 
lations.  but  in  the  body  stern  and  near-wake  regions  a  fifth- 
degree  polynomial  is  used  to  connect  the  two  displacement 
surfaces.  The  iteration  procedure  starts  with  a  potential  flow 
calculation  over  the  body  surface,  and  continues  until  the 
maximum  change  in  the  pressure  coefficient  on  the  displacement 
surface  is  less  than  a  specified  tolerance,  or  the  number  of 
iterations  is  exceeded.  The  results  converge  in  an  oscillatory 
manner,  and  so  the  final  pressure  coefficient  is  usually  taken  to 
be  the  average  of  the  previous  two  values.  The  original  version 
of  the  method  has  been  improved  in  two  ways  C5D.  Firstly,  a 
correction  was  made  to  the  predicted  velocity  profiles  to  take 
account  of  the  variation  of  static  pressure  across  the  boundary 
layer  and  to  obtain  the  correct  conditions  at  the  freestream 
boundary.  Secondly,  the  turbulent  mixing-length  formulation  was 
modified  to  include  thick  boundary- layer  effects.  Both  modifi¬ 
cations  are  aimed  at  giving  improved  predictions  in  the  propulsor 
region  at  the  stern  of  the  body. 


3.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 

In  this  Section,  brief  details  are  given  of  boundary- layer 
measurements  on  three  axisymmetric  bodies,  and  the  measurements 
are  compared  with  integral  and  finite-difference  method  predic¬ 
tions.  The  body  shapes  are  a  prolate  spheroid  with  conical  tail, 
which  was  investigated  by  Patel  et  al  CIOD.  and  two  in-house 
bodies:  a  torpedo-shaped  body  with  an  inflected  stern,  known  as 
the  propulsor  research  vehicle  (PRV)  C113  and  a  body  with  a 
streamline  nose  and  tail  cone  with  15-degree  semi-angle  C123. 
Figure  1  shows  the  geometry  and  notation  in  the  thick  boundary 
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layer  near  the  tail  of  a  typical  body,  and  the  body  shapes 
actually  used  are  plotted  as  Figure  2. 

The  modified  spheroid  was  investigated  experimentally  by 
Patel  et  al  CIOD  in  a  low  speed  wind  tunnel  at  Iowa  Institute  of 
Hydraulic  Research.  The  tunnel  speed  was  40  ft/sec  (12,2  m/s) 

corresponding  to  a  body-length  Reynolds  number  of  1.26x10^. 
Velocity  and  Reynolds  stress  measurements  were  conducted  usina  a 
hot-wire  anemometer,  with  the  boundary- layer  traverses  being 
normal  to  the  hull  enabling  direct  comparison  with  boundary-layer 
calculation  methods.  In  the  present  investigation,  predictions 
of  the  boundary- layer  integral  parameters  and  velocity  profiles 
were  obtained  using  both  integral  and  finite-difference  methods. 
Figures  3  and  4  show  comparisons  with  theory  of  the  measured 
integral  parameters  derived  by  integrating  the  velocity  measure¬ 
ments  across  the  boundary  layer.  Figure  5  shows  comparisons  of 
the  tangential  and  normal  velocity  components  as  predicted  by  the 
integral  and  finite-difference  methods  with  hot-wire  measure¬ 
ments.  The  integral  method  can  only  predict  the  tangential 
component.  The  finite-difference  predictions  here  do  not  take 
into  account  the  variation  of  static  pressure  across  the  boundary 
layer.  The  Integral  method  predictions  agree  well  with  the 
measurements  for  X/Ei<0.93.  At  X/I.-0.99.  the  finite-difference 
method  is  clearly  superior  to  the  integral  method  in  predicting 
the  tangential  velocity  component.  Figure  6  shows  similar 
comparisons  with  experiments  of  finite-difference  predictions, 
both  uncorrected  and  corrected  to  allow  for  static  pressure 
variation  across  the  boundary  layer.  The  improvement  between  the 
uncorrected  and  corrected  profiles  is  only  marginal  for  the 
tangential  velocity  component,  but  is  more  significant  for  the 
normal  velocity  component,  particularly  near  the  edge  of  the 
boundary  layer.  Figure  7  shows  measurements  and  finite-differ¬ 
ence  predictions  of  the  Reynolds  shear  stress.  The  agreement  is 
good  up  to  X/L=0.90  but  then  falls  off.  This  illustrates  the 
difficulty  of  using  a  simple  eddy-viscosity  model  of  turbulence 
in  a  thick,  axisymmetric  boundary  layer. 

The  propulsor  research  vehicle  has  a  nose  flat  and  both 
prediction  methods  were  modified  slightly  to  account  for  this. 
Velocity  measurements  were  carried  out  in  the  ARE  30-inch  Water 
Tunnel  at  Teddington  at  a  tunnel  speed  of  11  m/s  corresponding  to 
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a  body  length  Reynolds  number  of  3x10  .  The  body  diameter  was 
0.324m  and  the  tunnel  has  a  slotted  wall  working  section  0.76m  in 
diameter,  so  there  may  have  been  some  small  blockage  effect. 
Boundary-layer  velocity  components  were  measured  using  a  two-axis 
laser-Doppler  system.  Two  traverses  were  undertaken,  one  moving 
in  a  horizontal  direction,  where  the  axial  and  radial  components 
were  measured,  and  the  other  in  a  vertical  direction,  where  the 
axial  and  circumferential  components  were  measured.  The  tur¬ 
bulence  levels  were  also  estimated,  but  these  were  found  to  be 
relatively  high  because  of  the  slotted  wall. 

The  measurements  of  axial  and  radial  velocity  components  are 
compared  with  finite-difference  predictions  in  Figure  8.  The 
measured  profiles  are  along  radii  and  the  predictions  along 
normals  to  the  hull.  However  examination  of  the  predictions 
suggest  that  the  velocity  profiles  do  not  change  significantly  in 
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the  distance  spanned  by  the  measured  boundary  layer  and  so  no 
correction  was  applied.  The  finite-difference  method  is  shown  to 
aive  reasonable  predictions  of  both  the  axial  and  radial  velocity 
components . 

The  body  with  the  IS-degree  tail  cone  was  tested  in  the  No  5 
Wind  Tunnel  at  British  Maritime  Technology  Ltd.  at  a  body  length  ; 

Reynolds  number  of  6.5x10  .  The  main  objective  of  the  exper¬ 
iments  was  to  study  the  wake  of  a  body  with  fins,  but  some 

preliminary  measurements  were  carried  out  on  the  unappended  body.  i 

Tbe  axial  and  radial  velocity  components  were  measured  traversing 
the  boundary  layer  along  a  radius.  The  measurements  were  taken 
at  one  axial  and  seven  radial  locations  and  a  complete  circum- 
ferental  traverse  was  obtained  at  1°  intervals.  The  purpose  of 
this  experiment  was  to  ensure  that  any  departures  from  axial 
symmetry  due  to  the  effects  of  the  supporting  strut  and  wires 
were  minimal.  Axial  velocities  were  measured  traversing  in  a 
radial  direction.  These  are  compared  with  finite-difference 
predictions  in  Figure  9  and  the  agreement  is  generally  favour¬ 
able. 


4 .  CONCLUSIONS  AND  FURTHER  RESEARCH 

The  Wang  &  Huang  finite-difference  and  Myring  integral 
methods  for  predicting  the  turbulent  boundary  layer  on  axi- 
symmetric  bodies  have  been  evaluated  against  each  other  and 
compared  with  experiment.  The  methods  are  found  to  be  of  similar 
accuracy  over  all  but  the  extreme  tail  end  of  a  practical  body 
shape.  However  the  finite-difference  method  is  marginally 
superior  in  this  region  because  it  takes  some  account  of  static 
pressure  across  the  boundary  layer.  The  integral  method  can  only 
predict  the  tangential  velocity  component,  whereas  the  finite- 
difference  method  can  predict  the  tangential  and  normal  or  axial 
and  radial  components  and  the  Reynolds  shear  stress.  However  the 
integral  method  is  easier  to  use  and  produces  much  less  output 
than  the  finite-difference  method.  The  finite-difference  method 
is  recommended  for  the  maiority  of  future  computations,  but  some 
preliminary  calculations,  for  example,  to  investigate  the  effect 
of  varying  a  parameter,  are  probably  best  performed  using  the 
integral  method. 

There  are  several  more  sophisticated  axi symmetric  body 
prediction  methods  now  available  and  a  review  of  these  has 
recently  been  carried  out  by  Patel  &  Chen  C13D.  For  our  current 
body  shapes  they  appear  to  be  only  marginally  more  accurate  than 
the  finite-difference  method  but  require  significantly  more 
computational  effort.  For  practical  vehicle  configurations 
including  control  fins  and  propulsor,  a  three-dimensional 
Reynolds-averaged  Navier-Stokes  calculation  is  necessary. 

Research  into  the  relevant  numerical  methods  is  now  in  progress 
at  ARE. 


D  J  Atkins  (SSO) 
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FiG.l  BODY  GEOMETRY 


(a)  MODIFIED  SPHEROID 


(b)  PROPULSOR  RESEARCH  VEHICLE 


(c)  BODY  VWTH  15  DEC.  TAIL  CONE 


FIG.  2  BODY  SHAPES  TO  SCALE 
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FiG.3  COMPARISON  OF  PREDICTED  AND  MEASURED 
BOUNDARY  LAYER  PARAMETERS  FOR  MODIFIED  SPHEROID 
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FIGA  FURTHER  COMPARISON  OF  PREDICTED  AND  MEASURED 
BOUNDARY  LAYER  PARAMETERS  FOR  MODIFIED  SPHEROI 
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INTEGRAL  METHOD  PREDICTIONS 
FINITE-DIFFERENCE  METHOD  PREDICTIONS 


X  measurements  <U/UINF> 
measurements  <V/UINF) 


FIG.  5  COMPARISON  OF  PREDICTED  AND  MEASURED  TANGENTIAL 
AND  NORMAL  VELOCITIES  ON  MODIFIED  SPHEROID 
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FIG  6  COMPARISON  OF  FINITE-DIFFERENCE  PREDICTIONS  AND 
MEASUREMENTS  SHOWING  EFFECT  OF  STATIC  PRESSURE  CORRECTION 

(MODIFIED  SPHEROID) 
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FIG.  7  COMPARISON  OF  PREDICTIONS  AND  MEASUREMENTS  OF 
REYNOLDS  SHEAR  STRESS  ON  MODIFIED  SPHEROID 
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FIG.  8  COMPARISON  OF  PREDICTED  AND  MEASURED 
AXIAL  AND  RADIAL  VELOCITIES  ON  PROPULSOR  RESEARCH  VEHICLE 
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FIG.  9  COMPARISON  OF  PREDICTED  AND  MEASURED 
AXIAL  AND  RADIAL  VELOCITIES 
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